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ABSTRACT 

We have explored the relationship between hard X-ray (HXR) emissions and 
Doppler velocities caused by the chromospheric evaporation in two XI.6 class 
solar flares on 2014 September 10 and October 22, respectively. Both events 
display double ribbons and Interface Region Imaging Spectrograph {IRIS) slit is 
flxed on one of their ribbons from the flare onset. The explosive evaporations 
are detected in these two flares. The coronal line of Fe xxi 1354.09 A shows 
blue shifts, but chromospheric line of C 1 1354.29 A shows red shifts during the 
impulsive phase. The chromospheric evaporation tends to appear at the front 
of flare ribbon. Both Fe xxi and C i display their Doppler velocities with a 
‘increase-peak-decrease’ pattern which is well related to the ‘rising-maximum- 
decay’ phase of HXR emissions. Such anti-correlation between HXR emissions 
and Fe xxi Doppler shifts, and correlation with C i Doppler shifts indicate the 
electron-driven evaporation in these two flares. 

Subject headings: Sun: flares — Sun: UV radiation — Sun: X-rays, gamma rays — 
line: profiles — techniques: spectroscopic 
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Introduction 


Solar flares are drastic explosive phenomena on the Sun. They are able to release a 
huge amount of energy erg) in a typical time scale of tens of minutes. Based on 

the standard model, the flare energy is transferred from the magnetic energy. Magnetic 
reconnection is thought to be the primary energy release mechanism that heats the 
plasmas and accelerates the bi-directi onal nonthermal e l ectrons in the solar atmosphere . 


This is known as the C SHKP model (jCarmichael 


Kopp &: Pneuman 


1964 : 


Sturrock 


1966 


Hiravama 


1974 : 


19761) . These nonthermal electrons which guided by the reconnected 


magnetic held lines, not only travel to the interplanetary space but also precipitate into the 
lower corona and upper chromosphere, where they lose energy to produce radiation through 


Coulomb collisions with the denser medium. This 


model for the hard X-ray (HXR) emission fiBrown 


L. This 

las b 

Brown 

1971 


las been known as the ‘thick- 


Svrovatskii fc Shmeleva 


arget ’ 


197211 . 


Observations sho w that only a small 


(EUV) radiation flEmslie et al 


1978 


T action of the energy i s lost through e xtreme-ultraviolet 


Milligan et al. 


29141 : 


Milligan 


2915 1. The bulk of 


the energy heats the local chromospheric material rapidly up to a typical temperature 
of ~19 MK. Then the resulting overpressure can drive the mass flow upward along the 


loop at speeds of a few hundreds kilometers per secon d. The 


loops in a process c 

ailed ‘chromospheric evapora 

;ion’ 

Fisher et al. 

1985a 

: Liu et al. 

2996; 

Ning et al. 

2999 


2913 


Milligan 


flBrown 


rot p l asmas hll the 


1973 


Ning fc Cao 


Acton et al. 


2919 


flari ng 


1982 


Zhang fc Ji 


29151 ). resulting into soft X-ray (SXR) emission rising up. Substantia l 




2996; 

Ning 

et al. 

2999; 

Ning & Cao 

2919 

Ning 

29 

1; 

Ning & Cao 

2911 

Nitta et al. 

2912; 

Zhang & Ji 

2913 

, EU 

V fe.g.. 

Doschek et al. 

19§0 

Feldman et al. 

1989; 

Antonncci e1 

; al. 

I982I; 

Ding et al. 

1996; 

Teriaca et al. 

2996 

Milligan 

et al. 

2996a 

b; 

Millie 

i-an X Dennis 

2999 


Veronig et al. 


Chen fc Ding 


2919; 


29131 ) and radio flAschwanden fc Ben7Jll995 


Li fc Din 


2911 


Karlickv 


Doschek et a 


1998 


Ning et al 


2913 


Brosius 


299911 emissions. 
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Previous observations reveal that HXR emission tends to rise with the double footpoint 


sources along the flaring loop legs, e ventually merging them into a 


same position of the loop-top s ource flLiu et ah 


Nine fc Cao 


201 


1 


Ning 


2011 


Ning fc Cai 


3 


2006 


Ji et al 


20071. 


sing 


2008 


e source at the 


Ning et al. 


2009 


20111) . This is because the dense materials 


from the chromosphere rise upward along the loops, following the movement of the HXR 
emission targets. The flare shows the mergence speed of around 200 km s“^. Such large 


value of the mass evaporation is cl early demonstrated by the blue s. 


observations of the coronal lines flFeldman e 


2005 


Teriaca et al. 

2006; 

2010 


Chen & Ding 

2010; 

to 

o 

S — 

Young et al. 

2015) 


3rosius &: Holman 


Li fc Ding 


2011 


al. 


2007- 


980 


rifts a t the EUV spec tral 


Ding et al. 


1996 


Doschek et al. 


Milligan & Dennii 

3 ^ 


201f 


Berlicki et al. 


2009 


Brosius 


2013 


Veronig et al. 


Tian et al. 


20151) . Spectral images exhi 


at the outsides of the flare ribbons flCzavkowska 


r it that the blue shifts tend to appear 


1999 


Li fc Ding 


20041) . Evaporation 


materials with high temperature rise upward to disturb the coronal plasma, which results 
into the radio emission suddenly suppressed on the radio dynamic spectra, especially at 


Ning et al. 


200911 


lAschwanden & Benz 

1995; 

Karlicky 

1998 


From the observations, there are two types of chromospheric evaporation. Evaporation 
is said to proceed ‘gently’ when the chromosphere plasmas lose energy yia a combination of 
radiation and low-velocity hydrodynamic expansion. Emission lines formed at temperature 


characteristic of the at mosphere from 


all ap pear blue shifted flMilligan et al. 


chromosphere 


2nn6b 


Brosius 


tirou 


2009 


1 transition region t o corona 


Raftery et al. 


2009 


Li &: Ding 


20111 ). Evaporation is regarded to proceed ‘explosively’ when the chromosphere is unable 
to radiate energy at a sufficient rate and consequently expands at high velocities into the 
overlying flare loops. The oyerpressure of evaporated material also drives low-velocity 
downward motion into the underlying chromosphere, in a process known as ‘chromospheric 
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cond e nsation’ flWii 


20061 : 


(Wii 

ser el 

al. 

1994; 

Gzavkowska 

1999; 

Kamio et al. 

2005; 

Del Zanna et al. 


Teriaca et al. 


20061) ■ In this case, emission lines formed at temperature characteristic 


of the upper chromosphere and transition region all ap 


Dear rec 


shifted, while 


lotter 


lines from the corona appear : 

lue s 

lifted ( 

Fisher 

al. 

1985a 

b; 

Mil 

igan et al. 

2006a 

Del Zanna et al. 

20Qb; 

Brosins 

20Q9; 

Rafterv et al. 

200f 

; 

Li & Ding 

2011 

). Spectral 


observations show the red shifted velocity of ~20 to 40 km s“^, which is an order smaller 
than the blue shifted value (~200 km s“^). This is because the plasma density of the 
underlying lower chromosphere is much higher than that of the overlying corona. 


Up to now, there are two viewpoints about how to drive the evaporation in the 
literatures. One is the electron-driven, while another is the thermal conduction driven. The 
former emphases th at the non-thermal e n ergy of nonthermal electrons play an importa nt role 


in the evaporation flFisher et al. 


1985alJb : 


Milligan fc Dennis 


the la tter focuses on the thermal energy directly driven flFisher et al 


2009 


Tian et al. 


1985al: 


2014al). while 


Falewicz et al. 


2009|). In this paper, using the observations from Fermi Gamma-ray Burst Monitor (GBM), 


Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESST), Atmospheric Imaging 
Assembly (AIA) aboard Solar Dynamics Observatory (SDO) and Interface Region Imaging 
Spectrograph (IRIS), we explore the relationship between HXR emissions and Doppler 
velocities caused by evaporation during the impulsive phase in two solar flares in order to 
detect the observation evidences of the electron-driven chromospheric evaporation. 


2. Observations and Data Analysis 
2.1. Observations 


Two XI. 6 solar flares are selecte d to study in this paper. Firstly they are well covered 
by the IRIS flDe Pontieu et al.l 1201411 spectral observations. IRIS slit is fixed on the flare 
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ribbon during the impulsive phase, which gives us an opportunity to detect the whole 
history of Doppler velocities caused by evaporation. Secondly, HXR emissions are also well 


observed by Fermi flMeegan et al 


200911 or RHESSI flLin et ah 


20021) . One event takes place 


in NOAA AR 12158 on 2014 September 10. It starts at 17:21 UT and reaches its maximum 
at 17:45 UT on the GOES SXR light curves. Another occurs in NOAA AR 12192 on 2014 
October 22. It starts at 14:02 UT and peaks at 14:28 UT on the SXR emissions. 


Fig. [T] shows the SDO/AIA f Lemen et al. 2012) 131 A images (a, c) and /R/S'/SJ 
images (b, d) of these two flares. The contours on the AIA 131 A images r epresent the 


line-o f-sight magnetic fields from Helioseismic and Magnetic Imager (HMI) fjSchou et al. 


2 OI 2 I ) aboard SDO. The levels are set at 800 (purple) and -800 (orange) G, respectively. 


IRIS slit is fixed along the solar North-South direction on one ribbon of 2014 September 10 
flare. The cadence is 9.4 s. For the 2014 October 22 flare, IRIS slit is along the 45 degree to 
the North-South direction. IRIS detects the flare ribbon in the ‘raster’ mode. Each raster 
has eight steps (marked by the number on the SJI 1330 A image). Each step has a cadence 
of 16.4 s and a distance of ~2". Thus each raster has a duration of ~131 s. 

Both events display the double ribbons at SJI 1400 A or 1330 A images, as shown in 
Fig. [1] (b, d). The 2014 September 10 flare shows one short ribbon around the positive filed 
region, while another is long with a curved shape around the negative hied region. This 
long ribbon is dynamic and propagating toward the South-East direction, subsequently 
crosses the IRIS slit during the hare impulsive phase. Fig. [2] gives the time evolution of 
this ribbon on AIA 1600 A images. The arrow marks the ribb on propagation direction. 


This ribbon also exhibits strong Quasi-Periodic Pulsations (i.e.. 


Li fc Zhan 


2015 


Li et al. 


2 OI 5 I) . Similar to the 2014 September 10 hare, one ribbon of the 2014 October 22 hare is 
around the positive hied region, and the other ribbon at the negative hied region is detected 
by the IRIS slit during the hare impulsive phase. 
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2.2. Spectral fitting 


Fig. |3]shows the IRIS spectral profiles of three windows at ‘1343’ (a), ‘Fe xii ’ (b), and 
‘O I ’ (c) for the 2014 September 10 flare. IRIS has a spectral scale of ~25.6 mA/pixel for 
these three windows. The spectral data has been hrstly calibrated and processed with the 
routines of iris_orbitvar_corr_12.pro and iris_prep_despike.pro in the solar software (SSW) 


package. The hrst routine is us ed to correc 


the spacecraft orbital variation (ITian et ah 


; IRIS spectral image de 


2014b 


Cheng et ah 


brmation caused by 


20151 ). The second one is 


a generalized despiking tool for IRIS data. It could identify and remove the bad pixels 
through the iterative approach. At the flare onset (i.e., 17:28:43 UT), the spectral window 
at ‘O I ’ is characterized by many narrow, bright emission lines from neutral and singly 


the broad line of Fe xxi 1354.09 A (seen also., 


2014a . 


2015 


Young et al 


Li et al. 

2015; 

Polito et al. 

2015; 

Tian et al. 


20151 ). which is a typical coronal line to be used to detect the 


chromospheric evaporation. However, these blending chromospheric emission lines must be 
extracted before determining the Fe xxi intensity. According to the characteristics of IRIS 
spectral data, three steps are followed. 

Firstly, to determine the line centers and widths (FWHM: the full width at half 
maximum) at ‘O i ’ window, including blending lines and Fe xxi. There are seven blending 
lines of Fe xxi, including Fe ii 1353.07, 1354.06 A, Si ii 1352.69, 1353.78 A, C 1 1354.29 A, 
and two unidentihed lines at 1353.40, 1353.61 A, such as marked by the red vertical ticks 
(except for C i) in Fig. [3] (c). These lines are bright in the active regions while Fe xxi is 
quiet. Therefore, their line centers and widths can be detected from the single-Gaussian 


htting in the active regions. During the solar flare, their centers and widt 
in a r ange to do spectral htting based on the previous observations (e.g.. 


IS are constrained 


Curdt et al. 


2001 


20041 ). Taking the emission line of Si ii 1352.69 A for example, its center and width during 





































the flare are same as the values from the active regions but constrained, such as line center 
at 1352.69±0.102 A (ranging from 1352.588 to 1352.792 A), and the maximum width of 260 


mA, as 
studies 


isted in table □□ C i 


Doschek et ah 


1975 


las the strongest emission among these seven lines. Prev ious 


Cheng et ah 


1979 


Mason et ah 


1986 


Innes et ah 


2003a 


M) 


have identified its rest wavelength at 1354.29 A and its narrow width, which is assumed with 
maximum value of 130 mA for the flare spectral fitting in this paper. Its line center is set 
at the range of 1354. 29ifc0.26 A for the flare spectral fitting . Based on the fact that Fe xxi 


is a br oa d line (see., 


2003a 


U 


Tian et al 


Dose. 


2014al. 


rek et al 


2015 


1975 


Li et al 


Cheng et al 


1979 


Mason et al. 


19861: 


Innes et al. 


2015). Its line center is set as 1354.29±1.28 A, 


which almost cover the whole ‘O i ’ window, and its line width is assumed with a minimum 
value of 230 mA, as listed in Table [1] 


Secondly, to tie the blending line intensities from other similar isolated lines during the 
flare. Fig. |3] (c) shows coronal line of Fe xxi 1354.29 A blending with these seven emission 
lines for the 2014 September 10 flare. There is no way to determine their flare intensities 
only at ‘O I ’ window. However, IRIS has the spectra at other windows, i.e., ‘1343’, ‘Fe xii 
’. They also have the emission lines which behave similarly to the blending lines at ‘O i ’ 
window, such as H 2 1342.83 A at ‘1343’ window. Si ii 1350.13 A at ‘Fe xxi ’ window, and 
Fe II 1354.80 A at ‘O i ’ window. Their intensities can be used to tie the intensities of the 
six blending lines at ‘O i ’ window during the flare eruption, as listed in Table [1] This is 
because Si ii 1352.69, 1353.78 A at ‘O i ’ window have a similar behavior as Si ii 1350.13 A 
at ‘Fe XII ’ window, and Fe ii 1353.07, 1354.06 A at ‘O i ’ window have a similar behavior as 
Fe II 1354.80 A at ‘O i ’ window, and Unknown 1353.40, 1353.61 A have a similar behavior 
as H 2 1342.83 A at ‘1343’ window during the flare. These three emission lines are isolated 
and their intensities can be determined by a single-Gaussian fitting 


Thirdly, to determine the line parameters of Fe xxi and C i during the flare using the 
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multi-Gaussian fitting. Fig. |3] gives the six blending lines (except for C i) and three isolated 
emission lines marked with the red vertical ticks at three IRIS spectral windows. Fe xxi and 
C I are shown by the turquoise and magenta profiles. There are also another four isolated 
and bright lines marked with blue vertical ticks, such as Unknown 1348.34, 1348.60, and 
1349.65 A at ‘Fe xii ’ window, and Fe ii 1354.91 A at ‘O i ’ window. In total, these 15 
lines (i.e., H 2 1342.83 A, Unknown 1348.34 A, Unknown 1348.60 A, Unknown 1349.65 A, 
Si II 1350.13 A, Si ii 1352.69 A, Fe ii 1353.07 A, Unknown 1353.40 A, Unknown 1353.61 
A, Si II 1353.78 A, Fe ii 1354.06 A, Fe xxi 1354.09 A, C 1 1354.29 A, Fe ii 1354.80 A, Fe ii 
1354.91 A) superimposed on a linear background are used to do the multi-Gaussian htting 
at three IRIS spectral windows simultaneously, such as ‘1343’, ‘Fe xii ’ and ‘O i ’ windows. 
In our fitting method, both Fe xxi and G i have free intensities, almost free line centers and 
widths, as listed in Table [H The six blending lines of Fe xxi have constrained positions and 
widths and tied intensities. The other seven isolated lines have constrained positions and 
widths, but free intensities. Fig. |3] shows one example of such multi-Gaussian htting. The 
black prohles are the observational spectra at the positions of about 64.7" (marked by short 
orange line) on the slit for the 2014 September 10 hare. The brown prohles are 15 htting 
lines and the green line is the background. In this case, the intensities, centers and widths 
of 15 lines can be measured by the multi-Gaussian htting at any time and at any positions 
along the IRIS slit. Fig. [3] also shows that there are still some other unknown emission 
lines at these three spectral windows, such as 1342.09, 1344.08, 1348.03, 1350.75, 1352.02, 
1355.64 A. They are located at the edges of the spectral windows and do not ahect the 
spectral htting to detect Fe xxi and G i. Therefore they are not used to spectral htting in 
this paper. 

Fig. m gives four examples of the multi-Gaussian htting for the hares on 2014 September 
10 (top) and October 22 (bottom), respectively. The short vertical lines represent the 
rest wavelengths of Fe xxi (turquoise) and G i (magenta). It is well known that Fe xxi is 


10 


a hot coronal line with a formation temperatnre of abont 11 MK (log T 7.05), which 
resnlts into Fe xxi absent on the qniet Snn. Therefore, the rest wavelength of Fe xxi is not 


determined from the qniet regions. Recent studies from IRIS obse rvations show tha 


t Fe XXI 


Yonne; et al. 

2015; 

Graham & Ganzzi 

2015 


Tian et ah 


Polito et al. 


2014al. 


2015 : 


20151 ). In this paper, we use their 
average value of 1354.09 A as Fe xxi rest wavelength. C i is a typical chromosph eric line 


with a formation temperature of around 10^ K (log T 4.0) flHuang et al 


20141) . and its 


rest wavelength can be determined from the emissions at the quiet regions, as the dashed 
prohle shown in Fig. H] (b), which plots the spectral prohle average with 10 pixels around 
the position marked by the short black line. The spectra between 1353.66 A and 1354.68 A 
around C i are shown. In this paper, C i has a rest wavelength around 1354.29 A. 


2.3. Fitting parameters 

Using the method mentioned above, the intensities, line centers and widths of Fe xxi 
and C I can be determined from the multi-Gaussian htting along the IRIS slit. Fig. [5] shows 
the space-time diagrams of the intensities (a, c) and Doppler velocities (b, d) of Fe xxi and 
C I from 17:12 to 17:58 UT for the 2014 September 10 flare. There are strong EUV line 
emissions at about 60" — 80" and 35" — 40" along the slit. These two regions correspond to 
two propagation fronts (around 120 arcsec and 100 arcsec along the slit at 17:26 UT) of the 
curved ribbon in Fig. [2] (b), and the north region is the flare ribbon marked by the arrow. 
The spectral prohles at two positions of 64.7" and 60.6" on the slit but three different time 
are given in Figs. [3]and|l] (a, b). The 2014 September 10 flare ribbon starts to cross IRIS 
slit at about 17:25 UT although there are weak emissions from 17:12 UT. The flare ribbon 
front displays a narrow point at the beginning, then expands rapidly to a wide range of 
~25" after 17:40 UT along the slit. As mentioned before, this process suggests the flare 
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ribbon propagation across the IRIS slit. 


Doppler velocity is detected by the htting line center separation from the rest 
wavelength, whatever Fe xxi or C i. Fig. [5] (b) shows that the flare ribbon starts bine shifts, 
then red shifts at coronal line of Fe xxi, which is dne to the chromospheric evaporation 
at the beginning of the flare, then the hot materials fall back to the chromosphere along 


the flare 


OOP after cooling. This is consisten t with the standard flare model and previons 


hndings fjCzavkowska 


1999 


Li &: Ding 


20041) ■ and the evaporation appears at the onter 


edge of the flare ribbon. On the other hand, the evaporation tends to appear at the front 
of flare ribbon. The evaporation speed can reach ~230 km s“^, while the falling speed is 
~25 km s“^ in the corona. Fig. [5] (b) shows that the evaporation ronghly has a time scale 
of more than 10 minntes. Fig. [5] (d) gives the space-time diagram of C i Doppler velocities. 
Different from coronal line of Fe xxi, the chromospheric line of C i exhibits red shifts all 
dnring the flare, indicating that the 2014 September 10 flare is the explosive evaporation. 
There is a velocity peak of C i at the same time as Fe xxi blue shift. The velocity can reach 
~27 km s“^. After blue shift, Fe xxi displays the similar value of red shift as C i. 


Same as Fig. [5l Fig. [6] gives the space-time diagrams along the IRIS slit of the 
intensities and Doppler velocities of Fe xxi and C ifor the 2014 October 22 flare. As noted 
earlier, this event is observed in ‘raster’ mode with 8 steps. Thus, we can get 8 space-time 
diagrams at 8 step positions, respectively. Fig. [6] just shows the space-time diagram at the 
second step. In this case, the slit cadence is ~131 (16.4x8) s. The spectral profiles at two 
positions of about 59.2" and 47.4" on the slit are shown in Fig. 0] (bottom). Similar to the 
2014 September 10 flare, this event is the explosive evaporation. The coronal line of Fe xxi 
displays blue shift at the ribbon onset, then red shift, while the chromospheric line of C i 
exhibits red shift during the whole flare. The evaporation speed can reach ~145 km s“^, 
and the falling speed is around 20 km s“^. The evaporation timescale is roughly estimated 
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more than 10 minutes. 


Results 


In order to study the relationship between the HXR emission and the evaporation 
speed, Fig. [ 7 ] plots the X-ray light curves and the time evolution of Doppler velocities at Fe 
XXI and C i for both events. The 2014 September 10 flare is well detected by Fermi, but 
missed by RHESSI, which detects the 2014 October 22 flare. Fig. [ 7 ] (a) shows the GOES 
1 .0—8.0 A flux (black dashed lines) and Fermi/GBM light curves at 5 energy channels, such 
as 4.6-12.0 keV, 12.0-27.3 keV, 27.3-50.9 keV, 50.9-102.3 keV, and 102.3-296.4 keV. 
They are detected by the n2 detector, whose direction angle to the Sun is stable (~60°) 
before 17:45 UT, then its angle changes to the Sun to produce an X-ray peak, which is 
not real. There is a data gap after 17:54 UT. The time resolution of Fermi is 0.256 s, but 
becomes a higher value (0.064 s) automatically in the flare state. We rebin all the data 
into an uniform resolution of 0.256 s here, as shown in Fig. [ 7 ] (a). Fig. [ 7 ] (c, e) gives the 
time evolution of Fe xxi and C i Doppler velocities at two positions along the slit, i.e., at 
slit positions of 64.7" (orange) and 60.6" (purple) in Fig. [3]and|l] (upper). As predicted 
by the explosive evaporation, the coronal line of Fe xxi exhibits the blue shifts, while the 
chromospheric line of C i displays the red shifts at the same interval. Fe xxi increases its 
blue shifts rapidly to the maximum, then gradually and monotonically decreases to zero, 
continuously turns towards the red shifts. There are two possible explanations for these 
red shifts. Firstly, they maybe due to the hot material falling to the chromosphere along 
the flare loop after cooling. Secondly, they could be the signatures o 


seen in many imaging observations of flare arcades 


Sui fc Holman 


2003 


Li fc Gan 


2005 


2006 


Ji et al 




2006 


signal 

ures c 

loop contracting as 

Wang 

1992; 

Ambastha et ah 

1993 

, 2007: 

Zhou 

& Ji 

2009; 

Liu et al. 


2013 


Ning 


2013 


Yan et al. 


2013 


; IZhou et al. 


2013 


Kushwaha et al. 


2015 


Wang fc Liu 
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2015h . Similar to Fe xxi , C i increases its red shifts firstly to the maximum and then 
decreases to a stable state of about 24 km s There are two different physical scenarios 
to explain C i red shifts. The explosive peak of C i red shifts could be due to chromospheric 
condensation corresponding the evaporation detected as Fe xxi blue shifts at the same 
intervals, while the red shifts of C i on the decay phase could be due to the material falling 
back to chromosphere or the loop contracting, which results into the red shifts, not explosive 
but stable. Fe xxi and C i show explosive peak at the Doppler velocities, and the peak 
values can reach about -200 km s“^ and 27 km s“^, respectively. Although the peak time is 
different from the positions on the slit, the Doppler velocities exhibit the similar explosive 
peak. This is because the flare ribbon expands with time. The dashed lines in Fig. [7] (c) 
give the three times of the standard deviation (3cr) of the Doppler velocities from the quiet 
intervals (black prohles). The pluses (‘+’) mark the points where the speed values above 3cr 
and corresponding to the HXR peaks. Here, the evaporation time scale can be estimated 


from Fe xxi blue shifts. It is about 10 minutes for t 


le 2014 September 10 flare, which 


is con sistent with recent findings for the same flare (IGraham fc Canzzi 


2015 


Tianetah 


2015|). The pluses in Fig. [7] (e) mark the same points as that in Fig. [7] (c). After the blue 
shifts, Fe xxi becomes red shifts with a value of ~24 km s“^. Meanwhile, C i red shift 
velocity has the similar value (~24 km s“^), indicating the material falling downward with 
a similar speed from the corona. The error bars of the multi-Gaussian htting are displayed 
every 20-point with 2-5 uncertainties in Fig. [7] (c, e). The orange and purple colors are 
for two different positions on the slit, respectively. On the flare ribbon, the enhancement 
emissions result into the fitting speed errors decreasing to about 2 km s“^ (i.e., 5 = ~2 km 




. Same as Fig. [7] (a, c, and e). Fig. [7] (b, d, and f) shows the light curves of X-ray, Fe xxi 


, and C I Doppler velocities for 2014 October 22 flare, which is well detected by RHESSI. 
As mentioned before, this event is detected by IRIS in ‘raster’ mode, and Fig. [7] (d and 
f) shows the Doppler velocities at the position of the second step in each raster. Same 
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as Fig. [6l the time resolution is as low as 131 s. And the evaporation time scale of ~11 
minutes is estimated from Fe xxi blue shifts. The error bars of the multi-Gaussian htting 
are shown by every 2-point for each Doppler velocity. 

Fig. [ 7 ] (a) shows that there are three impulsive HXR peaks (> 27.3 keV) as marked by 
‘1’, ‘2’, and ‘3’. It is clear that the last two HXR peaks (‘2’ and ‘3’) are well correlated with 
the Doppler velocity peaks at two distinct positions, as shown in Fig. [ 7 ] (c, e), whatever Fe 
XXI or C I . In other words, both Fe xxi and C i exhibit a ‘increase-peak-decrease’ pattern 
of their Doppler velocities, well correlated with the ‘rising-maximum-decay’ phase at HXR 
emission. Considering the velocity direction, HXR light curves are anti-correlated with Fe 
XXI Doppler velocities, while correlated with C i Doppler velocities. This situation is well 
seen in the 2014 October 22 flare as well. There are two HXR peaks at the channel of 
25—50 keV, as marked by ‘1’ and ‘2’ in Fig. [ 7 ] (b). They are well correlated with the Fe xxi 
and C I Doppler velocity peaks at two distinct positions, as shown in Fig. [ 7 ] (d, f). 

Fig. [H] plots HXR peaks at 27.3—50.9 (or 25—50) keV dependence on the Fe xxi and C 
I Doppler velocities for the 2014 September 10 and 2014 October 22 flares, respectively. As 
expected from the electron-driven evaporation model, we hnd an anti-correlation between 
HXR emissions and coronal line (Fe xxi) evaporation velocities, while correlation between 
HXR emissions and chromospheric line (C i) condensation velocities. The correlation 
coefficient above 0.7 indicates that the non-thermal electrons cause the HXR emissions and 
drive the explosive evaporation simultaneously after precipitating in the chromosphere. 
Fig. m (c, d) shows only 4 points used for the correlation of each HXR peak due to the low 
time resolution of IRIS raster for the 2014 October 22 flare. For the peak ‘1’, the HXR 
emission at 25—50 keV starts an intensity as small as about 20 counts s“^, then becomes 
more than one order bigger after the maximum. That results into a single point at about 20 
counts s“^ of HXR emission in Fig. [8] (c, d). The correlation coefficient will become larger 
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if this point is omitted. 

Fig. Elplots the two HXR peaks at higher energy channels (such as 50.9—102.3 keV, and 
102.3—296.4 keV) dependence on Fe xxi and C i Doppler velocities for the 2014 September 
10 flare. The similar anti-correlation between HXR emissions and coronal line (Fe xxi) 
Doppler shifts, while correlation between HXR emissions and chromospheric line (C i) 
Doppler shifts are found during the same interval. The higher correlation coefficients (> 
0.85) further confirm our results. 


Discussions 


Using IRIS spectral observations on the flare ribbon and HXR observations from 
Fermi or RHESSI, we investigate the relationship between HXR emissions and Doppler 
velocities during the explosive evaporations in two X-class solar flares on 2014 September 
10 and October 22. Using the multi-Gaussian fitting, Doppler velocities of Fe xxi and C i 
are detected from IRIS spectral observations. At certain position on the slit, Fe xxi and 
C I display their Doppler velocities with a ‘increase-peak-decrease’ pattern, which is well 


findings flFisher et al. 


1985a 


b 


-decav’ phase of HXR emissions. Consistent with 

previous 

Milliean & Dennis 

2009; 

Brosius 

2013; 

Tian et al. 

2014a 

), 


we hnd a high anti-correlation between HXR emissions and coronal line Doppler shifts of Fe 
XXI , and a high correlation between HXR emissions and chromospheric line Doppler shifts 


of C I, indicating the electron beam-driven explosive ev aporation in 


The similar results are also found in the recent paper by 


Tian et al. 


hese two solar flares. 


fj2015l) . who get a high 


correlation between Fe xxi blue shifts and derivative of GOES SXR for the 2014 September 
6 and 10 flares. 


As noted earlier. Fig. [7] plots the Fe xxi and C i Doppler velocities at two distinct 
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positions respectively, which results into a good correlation with HXR peaks simultaneously. 
The peak time of Doppler velocities would be changed for various positions on the slit. 

In other words, Doppler velocities on the other positions are not well correlated with the 
HXR peaks at the same intervals. However, the time profiles of Doppler velocities at 
any positions exhibit the similar shape as shown in Fig. [7l We can still obtain the high 
correlations between Doppler velocities and HXR peaks after shifting an interval of Doppler 
velocity peaks. These intervals are different for the various positions along the slit, and it 
is an open question that a temporal correlation does not exist at these positions along the 
slit. Fig. [7]just gives the Doppler velocities at two special positions on the slit. They don’t 
need to shift with time to do the correlation in Figs. |8]and|9l For example, the first peak of 
the purple curve in Fig. [7] (e) is at 17:29 UT, and it could be correspond well with the peak 
‘2’ in Fermi data if it is shifted an interval of about 1 minutes. While th e second peak at 


17:33 UT are corresponding well with the HXR peak ‘3’ without shifting. 


Tian et al.l ( 2015 ) 


reported a time delay (~0.5—2.0 minutes) for the larger correlation between the Fe xxiblue 
shifts and SXR derivative on 2014 September 10 flare. This delay should be caused by the 
spectral prohles at different positions along the IRIS slit. When the flare ribbon propagates 
toward the South-West and expands with time, the peak of Fe xxi blue shifts changes with 
the slit positions, as well as with the time (seen in Figs. |2]and|5]). In this paper, we fit 
all the spectra along the IRIS slit at all time to obtain the space-time images of Doppler 
velocities, and get the time evolution of Doppler velocities at any positions on the slit. On 
the other hand, both Fe xxi and C i exhibit red shifts on the flare decay phase, especially 
for the 2014 September 10 flare. As the cross (‘x’) marked in Fig [7] (c, e), Fe xxi has a 
velocity of ~24 km s“^, while C i has a small peak. The spectral profile at this time is given 
in Fig. m (b), which shows strong emissions at Fe xxi and C i. At this position, Fe xxi and 
C I have line centers of about 1354.19 and 1354.40 A, line widths of 524.33 mA and 76.96 
mA. They show red shifts from their rest wavelengths. These red shifts of Fe xxi and C i 
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could be caused by the material falling down or the loop contracting at the decay phase of 
the flares. 

The rest wavelength of chromospheric line C i is set as 1354.29 A, which is detected 


from the quie 


observations fjPolito et ah 


Sun (seen in Fig. HI (b)). This is con sistent with recent studies from IRIS 


2015 


Sadvkov et ah 


2015|). However, as noted earlier, Fe xxi is a 


hot coronal line, which is absent on the quiet Sun. Therefore, we can not determine the rest 
wavelength of Fe xxi from non-flaring spectrum. There are different line centers of Fe xxi in 
the literatures . The first valu e of 13 54.1 A for Fe xxi has been identihed from the spectra of 
solar flares by 


Doschek et al 


1354.06 A and 135 4. 12 A fe.g.. 


limes et al. 


2003a 


b 


(1197511. Then t 


Cheng et al. 


Wang et al. 


le rest wavelength of Fe xxi is though between 


1979 


Mason et al. 


1986 


Feldman et al. 


2000 


200311 . In this paper, the rest w avelength o 


to 1354.09 A, anc 

this value is similar as 

Young et al. 

2015; 

Graham & Cauzzi 

2015 


hat in recent studies (iTian et al. 


Polito et al. 


2015 


Sadvkov et al 


Fe XXI is set 


2014a . 


2015 : 


2015) about 


Fe XXI from IRIS spectral data. Considering the broad range of Fe xxi, the rest wavelength 
has an uncertainty of about ±0.03 A, corresponding the Doppler velocity of about ±6.6 
km s“^. The rest wavelength of Fe xxi is probably taken from the decay phase of the flare 
due to it’s stable emission, as the example of the 2014 September 10 flare. However, the 
wavelength value is ~1354.19 A (i.e., the cross (‘x’) in Fig. [7] (c, d)) at the decay phase of 
this flare, which is much bigger than the rest wavelength in previous studies. It is a fact 
that the coronal line of Fe xxi is absent in the non-flare regions where the temperature is 
not hot enough. Therefore, the intensities and Doppler velocities of Fe xxi in Figs. E] (a, b) 
and [6] (a, b) outside the flare regions are invalid, they are fitting noises from the observation 
data. Same as the chromospheric line of C i, its Doppler velocities (Figs. |5] (d) and [6] (d)) 
outside the flare regions should also be fitting noises from the observation data. 


The multi-Gaussian fitting is used to determine the Doppler velocities of Fe xxi and 
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C I in this paper. There are three sources of uncertainties to their Doppler shifts. Firstly, 
as shown in Fig. [TJ the htting errors are very large in the non-flare regions and become 
much less in the flare regions, no matter Fe xxi or C i . This error is mathematic from 
the htting method. Secondly, the red wing enhancement or redshifted comp onent of the 


cool 

ines fi.e., F 

e II , Si II , and G i ) 

2015; 

Tian et al. 

2015 

). There is no 


Fe XXI at present time. Using tied intensities, line centers and widths could be a better 
approach as normally lines from the same ion or similar lines should have similar behaviors. 
Therefore, we constrained and tied these lines with the emission lines in other windows 
(seen in Table [1]) to eliminate the inhuence of these cool lines. Thirdly, the asymmetries 
of C I line should also affect the derived red shift of C i. Because the Gaussian htting of 


such asymmetrical lines tends to underestimate the v elocity of chromosp 


especially for the optically thick lines, i.e.. Mg ii line flGraham fc Cauzzi 


flDing et ah 


leric condensation. 


2015 '). and Ha line 


19951) . A goo d way to determ i ne the Doppler shi: 


the bisector method (i.e.. 


Ding et ah 


1995 


Graham fc Gauzzi 


s fro m the asymmetric line is 


20151) . which could accurately 


determine the Doppler shift, i.e., especially the red shift for the asymmetrical G i in our 
case. However, not the bisector, but the Gaussian htting is used in this paper. Because G i 
is one of the multi-Gaussian htting in our code, while the bisector method is a better way 
for the isolated line. 

In this paper, we hnd the blue shifts of Fe xxi quickly increases from z ero to the 
maxim um (more than 200 km s“^) before decreasing, which is diherent from 


Polito et ah 


(120151 ) hnding that the evaporation speed shows a monotonic decrease once appearing at 


the hare ribbon. This would be because we take the w hole history o 


velocities during the impulsive phase of the hare, while 


Polito et al. 


the Fe xxi Doppler 


(120151 ) only show the 


Doppler velocities after its maximum. On the other hand, their data is in ‘raster’ mode, 
rather than ‘sit and stare’ mode. The time resolution is not high. The 2014 October 22 hare 
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in this paper is in ‘raster’ mode, and has a lower time resolution than the 2014 September 
10 flare. But we can also detect the Dop pler velocity increase befo r e its maximum. T he flare 


on 2014 September 10 is also studied by 


Graham fc Cauzzi ( 2015 ): 


TianetaL 


fl2015h . Both 


of them found very clearly monotonic decrease of t he Fe xxi blue shifts, but no increase 


before the peak, as shown in Fig. 13 in 


Tian et ah 


(120151 ). They also found the evaporation 


within about 9 minutes, from 17:32 UT to 17:41 UT, which is similar to our finding of ~10 
minutes evaporation. From the orange curves in Fig. [7] (c), the start time of evaporation is 
about 17:26 UT (velocity above 3cr), and the end time is about 17 :36 UT (velocity equally 


to zero), while the peak time is about 17:28 UT. 


Tian et al, 


(120151) showed the monotonic 


decrease of Fe xxi blue shifts for the 2014 September 10 flare after the maximum at 17:32 
UT. Their velocity curve is at position around 117.8"on the slit. However, we give all the 
histories of the flare at the position of ~64. 7"along the slit in Fig. [7] (c), i.e., from the 


flare onset at 17:21 UT. On the other hand. 


Tian et al 


(120151) found the same evaporation 


pattern of ‘increase-peak-decrease’ for 2014 September 6 flare when their fitting results 
cover the whole histories of the flare. In fact, previous studies had been reported the 
similar evolution of Doppler velocities from hot coronal lines (i.e., Fe xii, Fe xix , Fe xxi 




so on) in t 

le solar flares, such as t’ 

e evolution o ’ 

‘increase-peak 

-decrease’ patt 

ern 

Wane et al. 

2003; 

Kamio et al. 

2005 

Brosius 

2009; 

Rafterv et al. 

2009; 

Li et al. 

2014 


Tian et al 


20151 ). Fig. [7| (e) also shows that the chromospheric C i line seems to display 


several small peaks after its maximum. They could be related to the HXR emissions at the 
decay phase or the materials with various speeds falling back to chromosphere. Meanwhile, 
they probably are the spectral signatures of the various loop contracting. 
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Table 1: The parameters of 15 emission lines at three IRIS spectral windows. 


IRIS window 

Ion 

Wavelength (A) 

Width (mA) 

Intensity tied to 


Si II 

1352.69±0.102 

< 260 

Si II 1350.13 


Fe II 

1353.07±0.051 

< 88 

Fe II 1354.80 


Unknown 

1353.40±0.061 

< 102 

H 2 1342.83 


Unknown 

1353.61±0.061 

< 102 

H 2 1342.83 

‘0 i’ 

Si II 

1353.78±0.102 

< 260 

Si II 1350.13 


Fe II 

1354.06±0.051 

< 88 

Fe II 1354.80 


Fe XXI 

1354.09±1.28 

> 230 



C I 

1354.29T0.26 

< 130 



Fe II 

1354.80±0.051 

< 88 



Fe II 

1354.91±0.061 

< 102 



Si II 

1350.13±0.102 

< 260 



Unknown 

1348.34±0.067 

< 102 


‘Fe XII ’ 

Unknown 

1348.60T0.067 

< 102 



Unknown 

1349.65±0.051 

< 77 


T343’ 

H2 

1342.83±0.061 

< 102 
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Fig. 1.— Top: AIA 131 A (a) and SJ 1400 A (b) images at 17:32 UT on 2014 September 10. 
Bottom: AIA 131 A (c) and SJ 1330 A (d) images at 14:12 UT on 2014 October 22. The 
bine and white lines represent the IRIS slit positions. The contours indicate the magnetic 
fields from HMI at the levels of 800 (purple) and -800 (orange) G. The red boxes mark the 
IRIS S JI regions. Two short blue lines on the IRIS slit in panels (c, d) mark the flare ribbon 
shown in Fig. [6l 
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Fig. 2.— Time sequences of AIA 1600 A images for the 2014 September 10 flare. The vertical 
line represents the IRIS slit positions, and the arrows mark the ribbon propagation. 
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Fig. 3.— Three IRIS spectra windows ((a) for ‘1343’, (b) for ‘Fe xii and (c) for ‘O i ’) at 
17:28:43 UT for 2014 September 10 flare. The black prohles are detected at ~64.7" along 
the slit positions (marked by the short horizontal line). The brown prohles represent the 
mnlti-Ganssian htting. The turqnoise prohle is Fe xxi , the magenta prohle is C i , and 
the green is the background. The other 13 emission lines used in this paper are labeled by 
vertical ticks. 
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Fig. 4.— Similar as Fig. [3] (c), IRIS flare spectra at ‘O i ’ window at four time on 2014 
September 10 and October 22, respectively. The black profiles are detected at various po¬ 
sitions on the slit marked by the short horizontal lines. The brown prohle represents the 
spectral htting. The turquoise prohle is Fe xxi , the magenta prohle is C i, and the green 
is the background. The short vertical lines mark the rest wavelengths of Fe xxi (turquoise) 
and C I (magenta), respectively. The dashed black prohle is the spectra from the non-haring 


region. 
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Fig. 5.— The space-time diagrams of intensities and Doppler velocities of Fe xxi and C ifor 
the 2014 September 10 flare. Y-axis is the distance along the whole IRIS slit. The spectral 
prohles at two positions on the slit marked by short horizontal lines and three times marked 
by vertical lines are shown in Figs. [3]and|l] (a, b) 
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Fig. 6.— Same as Fig.|5l but for the 2014 October 22 flare. Y-axis gives the distance between 


two short blue lines along the IRIS slit in Fig. [T] (c, d). The spectral prohles at two positions 
on the slit marked by short horizontal lines and two times marked by vertical lines are given 
in Fig. m (c, d). 
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Fig. 7 .— Top: the light curves from Fermi for the 2014 September 10 flare (a) and from 
RHESSI for the 2014 October 22 flare (b), and from GOES (dashed black lines) for both 
events. Middle and Bottom: the detected Doppler velocities of Fe xxi and C i at two dis¬ 
tinct positions on the slits (seen text in details). The dashed lines are the three times of 
the standard deviation, the pluses (‘+’) mark the points with speed values above 3a and 
corresponding HXR peaks. The cross (‘x’) marks the point in the decay phases of the flare, 
and its spectral prohle is given in Fig. |4] (b). The error bars represent the 2-6 uncertainties 
from the multi-Gaussian htting in the middle and bottom panels. 
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Fig. 8.— Scatter plots of two HXR peaks (orange and pnrple) dependence on Doppler 
velocities of Fe xxi and C i for the 2014 September 10 (a, b) and October 22 flares (c, d). 
The correlation coefficients (cc) are given. 
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Fig. 9.— Same as Fig. IH but HXR emissions at 50.9—102.3 keV and 102.3—296.4 keV for 


2014 September 10 flare. 






